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Abstract 

We describe propagation of torsional elastic waves in cylindrical waveguide with 
wedge dislocation in the framework of geometric theory of defects. The defect 
changes the dispersion relation. For positive deficit angles, it increases the phase 
velocity, and decreases the group velocity of the wave. 


1 Introduction 

Ideal crystals are absent in nature, and most of their physical properties, such as plasticity, 
melting, growth, etc., are dehned by defects of the crystalline structure. Therefore, a study 
of defects is a topical scientihc question of importance for applications in the hrst place. 
At present, a fundamental theory of defects is absent in spite of the existence of dozens 
of monographs and thousands of articles. 

One of the most promising approaches to the theory of defects is based on Riemann- 
Cartan geometry, which involves nontrivial metric and torsion. In this approach, a crystal 
is considered as a continuous elastic medium with a spin structure. If the displacement 
vector held is a smooth function, then there are only elastic stresses corresponding to 
diffeomorphisms of the Euclidean space. If the displacement vector held has disconti¬ 
nuities, then we are saying that there are defects in the elastic structure. Defects in the 
elastic structure are called dislocations and lead to the appearance of nontrivial geometry. 
Precisely, they correspond to a nonzero torsion tensor, equal to the surface density of the 
Burgers vector. Defects in the spin structure are called disclinations. They correspond to 
nonzero curvature tensor, curvature tensor being the surface density of the Frank vector. 

The idea to relate torsion to dislocations appeared in the 1950s [1-4]. This approach 
is still being successfully developed (note reviews [5-11]), and is often called the gauge 
theory of dislocations. 

Some time ago we proposed the geometrical theory of defects [12-14] . Our approach 
is essentially diherent from others in two respects. Firstly, we do not have the displace¬ 
ment vector held and rotational vector held as independent variables because, in general, 
they are not continuous. Instead, the triad held and §0(3)-connection are considered as 
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independent variables. If defects are absent, then the triad and §0(3)-connection reduce 
to partial derivatives of the displacement and rotational angle vector helds. In this case 
the latter can be reconstructed. Secondly, the set of equilibrium equations is different. We 
proposed purely geometric set which coincides with that of Euclidean three dimensional 
gravity with torsion. The nonlinear elasticity equations and principal chiral SO(3) model 
for the spin structure enter the model through the elastic and Lorentz gauge conditions 
[14-16] which allow to reconstruct the displacement and rotational angle vector helds in 
the absence of dislocations in full agreement with classical models. 

The advantage of the geometric theory of defects is that it allows one to describe 
singele defects as well as their continuous distributions. 

In the present paper, we consider propagation of torsional elastic waves in cylindrical 
waveguide with wedge dislocation. It is a classical problem which was solved in the 
absence of defect long ago within the elasticity theory (see, for example, m)- We show 
that the wedge dislocation changes the dispersion relation. For positive dehcit angles 
the phase velocity increases and group velocity decreases. Negative dehcit angles have 
opposite consequences. 

2 Elastic waves in media with dislocations 

Let us introduce notation and remind necessary facts from diherential geometry. We 
consider elastic media as topologically trivial manifold which is diheomorphic to Euclidean 
space Denote Cartesian coordinate system by x\ i = 1,2,3. Without dislocations, 
the elastic media is described by the Euclidean metric 6ij. If dislocations are present then 
the metric becomes nontrivial: 

^ij ^ ^ ^iji h) ^ 1) 2, 3, 

where e^*(x) is the triad. The triad held dehnes the orthonormal basis of the tangent 
spaces at each point e, := 6^*9^, where is the inverse triad. We shall see that this 
basis is more convenient in applications. 

Assume that relative displacements for elastic deformations are much smaller then 
deformations caused by defects: 

< e/. (1) 

where x) are components of the displacement vector held with respect to orthonormal 
basis in the tangent space e^. Then in the hrst approximation, the elastic waves propagate 
in Riemannian space with nontrivial metric produced by dislocations. Here we neglect 
changes in the metric produced by elastic waves themselves. Therefore for elastic waves 
in media with defects we postulate the following wave equation 

Poh* — /xAm* — (A + jj)V^V jU^ = 0, (2) 

where A and /i are Lame coefficients, A := V*Vj is the covariant Laplace-Beltrami op¬ 
erator for the triad held e^®, and V, is the covariant derivative. The explicit form of the 
covariant derivative is 


S/iU^ := (3) 

where is the §0(3)-connection constructed for zero torsion. 

Note that the displacement vector held u® is not the total displacement vector held of 
points in the media with dislocations. In those regions where defects are absent, the total 
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displacement vector is u\^ + m* where u\^ is the dislocation part of the displacement vecor 
held which produces the triad := d^u\y. Note also that the smallness of relative elastic 
deformations Q is meaningful at those regions where displacements uj, are not dehned 
(they are not continuous functions if dislocations are present). 

The expression for the deformation tensor cjjj is needed to pose the boundary value 
problem for wave equation ([^. In the presence of defects the deformation tensor is dehned 
as follows 


^ij ■ T fi^i) ■ 


(4) 


Thus the existence of dislocations in media results in the introduction of the triad held 
e^*(x), metric gfj,u{x), and replacement of partial derivatives by covariant ones. 

Remember that lowering of Latin indices is performed using the Euclidean metric, 
Ui '■= SijU^, and it commutes with covariant diherentiation. 


2.1 Torsional waves in cylindrical waveguide with wedge dislo¬ 
cation 

We consider cylindrical waveguide of radius a with a wedge dislocation. Its axis is assumed 
to coincide with the axis of cylinder. We choose cylindrical coordinate system r, 93 , z, the 
axis coinciding with the cylinder axis as well. Let the wedge dislocation has dehcit angle 
27i9 (see FigQ If —1 < 0 < 0, then the wedge is cut out from the cylinder. Dislocation 



Figure 1: Wedge dislocation with the dehcit angle 2710. For negative and positive 6 *, the 
wedge is cut out or added, respectively. 

is absent for 9 = 0. For positive dehcit angle, 9 > 0, the wedge of media is added to the 
cylinder. 

The cylinder contains wedge dislocation, and it creates nontrivial metrics. The metric 
becomes noneuclidean na 

where a := 1 + 9, and we introduce dimensionless constant 

7 ;= -eh + ^/9%'^ + 1 + 9. ( 6 ) 
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The constant 

b ■= . ^ . 

2(1 - a) 

is defined by dimensionless Poisson ratio 

A 

2(A + /i)’ 

which characterize elastic properties of media. Thermodynamical argnments restrict pos¬ 
sible valnes of Poisson ratio to—1<(T<1/2 [T8] . 

The limit 

a ^ 1, 7^1 (7) 

corresponds to absence of the dislocation. 

Metric (|^ is rotationally and translationally invariant. Two dimensional metric on 
sections 2 ; = const is the metric of conical singnlarity which arises dne to the wedge 
dislocation. It is important that this metric is written in elastic gange and therefore is 
consistent with the elasticity theory [T5] . 

Metric for the wedge dislocation (|^ is defined for all deficit angles 9 > —1 and all 
valnes of radial coordinate 0 < r < a. It is degenerate on the dislocation axis r = 0. 
Metric (j^ coincide with the indnced metric in classical elasticity theory only for small 
deficit angles and near the bonndary of the cylinder r ~ a prams]. 

Metric for the wedge dislocation defines the triad field which we choose to be diagonal: 


Cr = ( - ) 


r\ 7-1 
a) 


/r\ 7—1 ar 

‘P = (-] — e^ = l 

\a/ 7 


( 8 ) 


where = r,ip,z and i = r,ip, z. In the absence of dislocation, 7 = 1, a = 1, it defines the 
nsnal orthonormal basis for tangent spaces in cylindrical coordinates. The inverse triad 

is 


'a\7-i 
e' r = I - 


= I - 


a\7-i 7 


- 6^ = 1 

ar 


(9) 


To find the explicit form of the wave operator, we need to compnte Christoffel’s sym¬ 
bols and components of §0(3)-connection. Straightforward calcnlations show that only 
for Christoffel’s symbols differ from zero: 


P ^ = 

J. 


-p V — -p — ± 

J- 

r 


7-1 


7 


( 10 ) 


r — 

J- ipip 


o 

a r 
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Triad ([^ defines also components of §0(3)-connection 

^ J = -d e -I- f Pe ^ 


Only two components differ from zero: 

0 _ r _ 

— a. 


( 11 ) 


Now we can compnte the Laplacian 


Aui = g^'^dl^Ui - - 2gP'"u^id^Uj- 

- gp’'ff,/{dpUi - Upi^Uj) + gp’-'upi^ojuAuj. ( 12 ) 
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Substitution of explicit expressions for Christoffel’s symbols (10) and SO(3)-connection 
(0 yields the following expressions: 


Auf = ( - ) 
-r / 


a\27-2 1 ^ / „ ^ /a\27-2 ^2 


-dr{rdrUf) + 

27-2 ^2 




-I-) 

.r / 


—Mr — 2 - 


a\ 


27-2 ^2 


2 , 


ar 


Am<3 = I - 
.r / 


a\27-2 1 


^dr{rdrU^ A 


a\ 27-2 ^2 


q ,2^2^ w '“‘^ + dzzUif 


(13) 




27-2 ry 2 




27-2 ^2 


_ 2 ^‘P^ri 


ar 


Auz = 1 - ) 

r / 


a\27-2 1 


-dr{rdrUz) + [-j 


a\27-2 ^2 


d,„,„Uz “ 1 “ dzzUz 




Corresponding wave equations ([^ have many solutions. The simplest one describes 
torsional waves for which only the angular component differs from zero and does not 
depend on the angular coordinate (p: 


Ur = 0 , 


M,; 




Uz = 0 . 


It is easy to see that torsional waves take place without media compression 

e := ViM* = 0 . 

For these oscillations, wave equations \3ur = 0 and Auz = 0 are identically satisfied, and 
we are left with one wave equation 


~ Vr/ -dr[rdrUY - ~ ^ 

We look for solution of this equation in the plain wave form 


N 27-2 ^2 


(14) 


= re {U e" 


where U{r) G M is the wave amplitude, A: G M is the wave vector, and a; G M is the wave 
frequency. Then wave equation (14) takes the form 


rdr {rdrU) + j U — y^f/ = 0, 


where 


;= _ Y 


(15) 


(16) 


Now we introduce new radial coordinate 


r = ar 'y, 


0 < r' < 1 . 


Then equation (16) reduces to the Bessel equation 

12 d U I dU X 2 /2 


r ■ 


dr''^ 


+ r'— + W^U-U = 0, 
dr 


(17) 
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where 


2^2 


A" := 


n a 


T 


A general solution to the Bessel equation has two integration constants. We require it to 
be hnite at r = 0. Then the amplitude takes the form 




(18) 


where A G M is the integration constant and Ji is the Bessel function of the first kind of 
order one (see, for example, IZDl). 

If dislocation is absent, 7 = 1 , then the solution takes the well known form (see, for 
example, irzi) 

U AJy(Kr) 

7^1 

To impose boundary conditions on the cylinder surface, we need explicit form of de¬ 
formation tensor (|^. Straightforward calculations yeild the following expressions: 


^rr — ^ ^ 

= (“) ^ auf), 


1 /a\7-i 


“ 2 U 


^ZZ dz'llz') 

Only two components 


= 2 
_ 1 
- 2 


1 /a\T'-i 


a\7-i 


7 7 

d’pUin -f" dmUr U(p 

ar r 


drUz + dzUf 


(19) 


^ 1 
1 \r 


-) —d^Uz + dzU^ 
r / ar 


a ^ ' 

Uip 

r 


differ from zero for torsional waves. We require the cylinder surface to be free. Then the 
elastic forces on the boundary must be zero. It yeilds the boundary condition 


o ^ ' 

Uj'ZLfn U0 

r . 


( 20 ) 


For solution (18), it takes the form 


AJ((A)-Ji(A) = 0, 

where prime denotes differentiation of the Bessel function with respect to its argument. 
There is equality 

j:(a) = MX) - 7 .(a). 

Then the boundary condition is equivalent to the equality which dehnes the dispersion 
relation 


7 \ 


( 21 ) 


where ^ is a root of the equation 


eJo(0 = 2Ji(0. 


( 22 ) 
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Bessel functions have the following asymptopics for large argument, ^ S> 1, ^ ^ z/, 


Jo(A) ^ 






Ji(A) 




Then equation (21) takes the form 


TT 


K - T = 2sin u - T 


TT 


for large argument. This equation has countable set of roots, each of them dehning the 
dispersion relation. 

The phase velocity of torsional waves v \ = uj/k has the form 




V = Cxy 1 + 

It is not difficult to calculate the group velocity 

doj 

^ dk V 


(23) 


(24) 


If dislocation is absent, equation (21) for ^ remains the same. Relation (24) between 
phase and group velocities preserves its form. Thus the presence of the wedge dislocation 
changes only dispersion relation (21). 

For positive dehcit angles 7 > 1, the phase velocity increases and the group velocity 
decreases as the consequence of equations (23) and (24). 

For small dehcit angles the hrst correction has the form 

1 - 2a 


7 1 + 6* 


2(1 


cr 


and 


u 




i + e 


1 - 2a 


2(1-a) 1 + ^ 


(25) 


It is linear in dehcit angle of the dislocation. 


3 Conclusion 

We showed that the presence of the wedge dislocation in cylindrical waveguide leads to 
changing of the dispersion relation for torsional elastic waves. The diherence increases 
as long as the dehcit angle increases. For positive dehcit angle, the phase velocity of the 
wave increases and group velocity decreases. For negative dehcit angle, the behavior is 
opposite. 

Solution of this problem within the classical elasticity encounters essential technical 
problems. Indeed, relative displacements are not small in the neighbourhood of the dis¬ 
location axis even for small dehcit angles, and the expression for the induced metric is 
much more complicated [191 IIS]- Solution of this problem within the geometric theory of 
defects as simpler. 

This work is supported by the Russian Science Foundation (project 14-11-00687) in 
Steklov Mathematical Institute. 
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